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Two-particle azimuthal (A^) and pseudorapidity (Ary) correlations using a trigger particle with 
large transverse momentum (pr) in d+Au, Cu+Cu and Au+Au collisions at y^sjvjv ~ 62.4 GeV and 
200 GeV from the STAR experiment at RHIC are presented. The near-side correlation is separated 
into a jet-like component, narrow in both A^ and A77, and the ridge, narrow in A0 but broad in 
Arj. Both components are studied as a function of collision centrality, and the jet-like correlation 
is studied as a function of the trigger and associated px- The behavior of the jet-like component 
is remarkably consistent for different collision systems, suggesting it is produced by fragmentation. 
The width of the jet-like correlation is found to increase with the system size. The ridge, previously 
observed in Au-|-Au collisions at ^sjvjv = 200 GeV, is also found in Cu-|-Cu collisions and in collisions 
at y/SNN ~ 62.4 GeV, but is found to be substantially smaller at ^sjviv = 62.4 GeV than at ^/sNM 
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— 200 GeV for the same average number of participants ((A'pa 
compared to models. 

PACS numbers: 25.75.-q,21.65.Qr,24.85.+p,25.75.Bh 



t)). Measurements of the ridge are 
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107 I. INTRODUCTION 

108 Jets are a useful probe of the hot, dense medium ere-"" 

109 ated in heavy-ion collisions at the Relativistic Heavy 
no Ion Collider (RHIC) at Brookhaven National Laboratory"^ 

111 (BNL). Jet quenching [l[ was first observed as the sup-"^ 

112 pression of inclusive hadron spectra at large transverse 

113 momenta (pt) in central Au+Au collisions with respect"'^ 

114 to p+p data scaled by number of binary nucleon-nucleon"'* 

115 collisions 0-13 • Properties of jets at RHIC have been"^ 

116 studied extensively using di-hadron correlations relative"** 

117 to a trigger particle with large transverse momentum 

118 [a. 

119 Systematic studies of associated particle distributions 

120 on the opposite side of the trigger particle revealed their 

121 significant modification in Au+Au relative to p+p and 

122 d+Au collisions at the top RHIC energy of ys/v7V = 200 

123 GeV. For low passociatod^ ^Yic amplitude of the away-side 

124 peak isgreater and the shape is modified in Au-|-Au colli- 

125 sions At intermediate pr (4 < p^r^'^^"' < 6 GeV/c,"' 

126 2 GeV/c < pl^'^o^ated ^ p^"gg<=^)^ the away-side corre-"' 

127 lation peak is strongly suppressed At higher Pt,„^ 

128 the away-side peak reappears without shape modifica-^^^ 

129 tion, but the away-side per trigger yield is smaller in^^^ 

130 Au-h Au collisions than in p+p and d-f- Au [ll| . 

131 The associated particle distribution on the near side ofiai 

132 the trigger particle, the subject of this paper, is also sig-i82 

133 nificantly modified in central Au+Au collisions. In p+Pis^ 

134 and d+Au collisions, there is a peak narrow in azimuthi84 

135 (A0) and pscudorapidity (Ayy) around the trigger par-i85 

136 tide, which we refer to as the jet-like correlation. ThiSi86 

137 peak is also present in Au+Au collisons, but an addi- 

138 tional structure which is narrow in azimuth but broad 

139 in pseudorapidity has been observed in central Au+Au 

140 collisions at ^/sWn = 200 GeV 0, This struc-^^^ 

141 turc, called the ridge, is independent of At; within the^^^ 

142 STAR acceptance, \Ar]\ < 2.0, within errors, and persists^^^g 

143 to high pT (p^"*^*^"' « 6 GeV/c, ^associated ^ 3 GeV/c).i<,„ 

144 While the spectrum of particles in the jet-like correla-m 

145 tion becomes flatter with increasing py'^*""^, the slope^gs 

146 of the spectrum of particles in the ridge is independentj(,3 

147 of py'^^'^' and closer to the inclusive spectrum than to^g^ 

148 that of the jet-like correlation. Recent studies of di-^gj 

149 hadron correlations at lower transverse momenta (p^'^^°' 

150 > 2.5 GcV/c, p|?>50"atcd > 20 McV/c) by the PHOBOSig, 

151 experiment show that the ridge is roughly independcntigg 

152 of A?; and extends over four units in A77 [isj . A sim-^g, 

153 ilar broad correlation in pseudorapidity is also evidentjoo 

154 in complementary studies of minijets using untriggeredsoi 

155 di-hadron correlations [l6l . [iTj . 202 

156 Several mechanisms for the production of the ridge203 

157 have been proposed since the first observation of this204 



new phenomenon. In one model [18| the ridge is pro- 
posed to be formed from gluon radiation emitted by a 
high-pT parton propagating in the medium with strong 
longitudinal flow. The momentum-kick model proposes 
that the ridge forms as a fast parton traveling through 
the medium loses energy through collisions with partons 
in the medium, causing those partons to be correlated 
in space with the fast parton jl9l - l2]| . Parton recombi- 
nation has been also proposed as a mechanism for the 
production of the ridge [23 - l23 | . Another model [2^, 
suggests that the ridge is not actually caused by a hard 
parton but is the product of radial flow and the surface 
biased emission of the trigger particle, causing an ap- 
parent correlation between particles from the bulk and 
high-pT trigger particles. 



Another class of models is based on the conversion 
of correlations in the initial state into momentum space 
through various flow effects. The model in (23] explains 
the ridge as arising from the spontaneous formation of ex- 
tended color fields in a longitudinally expanding medium 
due to the presence of plasma instabilities. Long-range 
pseudorapidity correlations formed in an initial state 
giasma combined with radial flow have been also dis- 
cussed as a mechanism for the ridge [28l - [30j . Recently, 
it has been suggested that triangular anisotropy in the 
initial collision geometry caused by event-by-event fluc- 
tuations can give rise to triangular flow, which leads to 
the ridge and contributes to the double peaked away-side 
observed in heavy- ion collisions at RHIC (3ll - [37| . 



In this paper we present measurements of the sys- 
tem size and collision energy dependence of near-side di- 
hadron correlations using data from Cu+Cu and Au+Au 
collisions at ^/s^n = 62.4 GeV and ^/JWn = 200 GeV 
measured by the STAR experiment at RHIC. In partic- 
ular, we investigate the centrality dependence of the jet- 
like correlations and the ridge and the transverse momen- 
tum dependence of the jet-like correlations. The prop- 
erties of jet-like correlations in heavy-ion collisions are 
compared to those from d+Au collisions and PYTHIA 
simulations to look for possible medium modifications 
and broadening of the near-side jet-like correlation, for 
example, due to gluon bremsstrahlung [isj . The new re- 
sults on the system size and energy dependence of the 
ridge yield presented in this article extend our knowl- 
edge of this phenomenon and, in combination with other 
measurements at RHIC, provide important quantitative 
input and constraints to model calculations. 
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II. EXPERIMENTAL SETUP AND DATA 
SAMPLE 



The results presented in this paper are based on data 
measured by the STAR experiment from d+Au colhsions 
at y^JWN = 200 GcV in 2003, Au+Au colhsions at y^JWW^^^ 
= 62.4 GeV and 200 GcV in 2004, and Cu+Cu collisions 
at y^IWN = 62.4 GeV and 200 GeV in 2005. The d+Au^" 
events were selected using a minimally biased (MB) trig- 
ger requiring at least one beam-rapidity neutron in the254 
Zero Degree Calorimeter (ZDC), located at 18 m from255 
the nominal interaction point in the Au beam direction, 256 
accepting 95±3% of the Au-(-Au hadronic cross section257 
[ssf . For Cu-|-Gu collisions, the MB trigger was based258 
on the combined signals from the Beam-Beam Gounters259 
(BBC) at forward rapidity (3.3 < \t]\ < 5.0) and a coin-26o 
cidence between the ZDCs. The MB trigger for Au-|-Au26i 
collisions at -y/sTvw = 62.4 GeV and 200 GeV was ob-262 
tained using a ZDC coincidence, a signal in both BBCs263 
and a minimum charged particle multiplicity in an array264 
of scintillator slats arranged in a barrel, the Central Trig- 
ger Barrel (CTB), to reject non-hadronic interactions. 

For Au+Au collisions at -y/s^vjv = 200 GcV, an addi- 
tional online trigger for central collisions was used. This 
trigger was based on the energy deposited in the ZDCs 
in combination with the multiplicity in the CTB. The 
central trigger sampled the most central 12% of the totabisB 
hadronic cross section. 266 

In order to achieve a more uniform detector accep-^**^ 
tance, only those events with the primary vertex position^''*' 
along the longitudinal beam direction (z) within 30 cm^''' 
of the center of the STAR detector were used for the^™ 
analysis. For the d-l-Au collisions this was expanded to^" 
l^l < 50 cm. The number of events after the vertex cut^^^ 
in individual data samples is summarized in Table HI 

The STAR Time Projection Chamber (TPC) ^ was 
used for tracking of charged particles. The collision cen- 
trality was determined from the uncorrected number of 
charged tracks at mid-rapidity (I77I < 0.5) in the TPC."" 
The charged tracks used for the centrality determina-"' 
tion had a three dimensional distance of closest approach"*" 
(DCA) to the primary vertex of less than 3 cm and at"^ 
least 10 fit points from the TPC. Each data sample was"" 
then divided into several centrality bins, and the fraction"' 
of the geometric cross section, the average number of par-^™ 

281 



TABLE I: Number of events after cuts (see text) in the 
samples analyzed. 
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System 


Centrality 


^SNN [GeV] 


No. of events [lO''] 


Cu-^Cu 


0-60% 


62.4 


24 


Au-I-Au 


0-80% 


62.4 


8 


d+Au 


0-95% 


200 


3 


Cu+Cu 


0-60% 


200 


38 


Au+Au 


0-80% 


200 


28 


Au+Au 


0-12% 


200 


17 



290 
291 
292 
293 



ticipating nucleons ((A^part)) and the average number of 
binary collisions ((A'coii)) were calculated using Glauber 
Monte-Carlo model calculations 1401. 



III. DATA ANALYSIS 

A. Correlation technique 

Tracks used in this analysis were required to have at 
least 15 fit points in the TPC, a DCA to the primary 
vertex of less than 1 cm, and a pseudorapidity \ri\ < 1.0. 
As in previous di-hadron correlation studies in the STAR 
experiment [H, |4l|, 14^ , a high-pr trigger particle was se- 
lected and the raw distribution of associated tracks rela- 
tive to that trigger in pseudorapidity (At;) and azimuth 
(A^) is studied. This distribution, d^A^raw/^A^dA?;, is 
normalized by the number of trigger particles, A'trigger, 
and corrected for the cfSciency and acceptance of associ- 
ated tracks e: 



dAcj) dAr] 



(A0,A77) 



1 



A^triggor dA(t)dAr] 
1 1 



b.-q) epair(A0, A-q) 



(1) 



The efficiency correction £assoc(0, is a correction for 
the TPC single charged track reconstruction efficiency 
and epair(A(^, Ar;) is a correction for track merging and 
finite TPC track-pair acceptance in A^ and A77 as de- 
scribed in detail below. The data presented in this paper 
are averaged between positive and negative At/) and At; 
regions and are reflected about A(f) =0 and At; =0 in the 
plots. 



B. Single charged track efficiency correction 

The single charged track reconstruction efficiency in 
the TPC is determined by simulating the detector re- 
sponse to a charged particle and embedding these signals 
into a real event. This hybrid event is analyzed using the 
same software as for the real events. The efficiency for 
detecting a single track as a function ofpr, ??, and central- 
ity is determined from the number of simulated particles 
which were successfully reconstructed. The single track 
efficiency is approximately constant for px > 2 GeV/c 
and ranges from around 75% for central Au+Au events 
to around 85% for peripheral Cu+Cu events as shown 
in Figure [TJ The efficiency for reconstructing a track in 
d+Au is 89%. The systematic uncertainty on the effi- 
ciency correction, 5%, is strongly correlated across cen- 
tralities and pt bins for each data set but not between 
data sets. In the correlations, each track pair is corrected 
for the efficiency for reconstructing the associated parti- 
cle. Since the correlations are normalized by the number 
of trigger particles, no correction for the efficiency of the 
trigger particle is necessary. 
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FIG. 1: (Color online) Parameterizations of the transverse 
momentum dependence of the reconstruction efficiency of^"" 
charged particles in the TPC in various collision systems, en-^"^ 
ergies and centrality bins for the track selection cuts used in^"^ 
this analysis. Note the zero suppression of the axes. 



294 C. Corrections for track merging and track 

295 crossing effects in the TPC '47 

348 

296 The reconstruction of charged tracks from TPC hits is^"' 

297 performed iteratively, with hits removed from the event 

298 once they are assigned to a track. If two tracks have small 

299 angular separation in both pscudorapidity and azimuth, 

300 they arc more difficult to reconstruct because distinct hits 

301 from each particle may not be resolved by the TPC. If two 

302 particles are close in momenta or have sufficiently high 

303 pt that their tracks are nearly straight, they may not be 

304 distinguished. This effect, called track merging, reduces 

305 the number of pairs observed at small opening angles and 

306 results in an artificial dip in the raw correlations centered^''' 

307 at (A<?i,Ar/) = (0,0). 

308 Figure [2] shows an example of a {A(f>, Arj) two-particle 

309 di-hadron correlation function in central Au+Au colli- 

310 sions at ^/snn = 200 GeV, along with the corresponding^^^ 

311 Acj) and A77 projections for 3 < < 6 GeV/c andj^^^ 

312 1.5 GeV/c < p|?>5°"atcd < pWggor g^^y angular sep-366 

313 arations, a clear dip in the raw correlations is visible and367 

314 must be corrected for in order to extract the yield of36s 

315 associated particles on the near side. 369 

316 Another similar effect is evident in the active TPC vol-370 

317 umc from high-pr tracks which cross. While few hits are37i 

318 lost in this case, one track may lose hits near the cross-372 

319 ing point and therefore be split into two shorter tracks. 373 

320 Shorter tracks are less likely to meet the track selection374 

321 criteria. Track merging and track crossing cause four375 

322 dips near (Acf), A77) = (0,0) but slightly displaced in A0.376 

323 The location and width of these dips in Acf) is dependent377 

324 on the relative helicities, h, and the px intervals of the378 

325 trigger and associated particles. The hclicity h is given379 



h = 



-qB 



(2) 



where q is the charge of the particle and B is the magnetic 
field. The dips for tracks of the same hclicity arc dom- 
inantly due to track merging and the dips for tracks of 
opposite helicities are dominantly due to track crossing. 
Figure [3] displays the correlation function from Figure [2] 
in four different hclicity combinations of trigger and asso- 
ciated particles showing the finer substructure of the dip 
on the near side. When the helicities of the trigger and 
associated particles are the same, the percentage of over- 
lapping hits is greater. Because higher pT tracks have a 
smaller curvature, it is more likely for two high pT tracks 
close in azimuth and pscudorapidity to be merged than 
lower Pt particles. However, track pairs are lost whether 
the pair is part of the combinatorial background or part 
of the signal. This effect means that the magnitude of the 
dip is greater in central collisions where the background 
is greater and decreases with increasing pT because of the 
decreasing background. 

One way to correct for the track merging effect is to 
remove pairs from mixed events that would have merged 
in real data. The environment in mixed events must be 
similar to real data in order for pair rejection to be ac- 
curately reproduced. The reference multiplicities of both 
events were required to be within 10 of each other to 
assure similar track density. In addition, mixed events 
were required to have vertices within 2 cm of each other 
along the beam axis in order to ensure similar geometric 
acceptance to avoid a different dip shape in A77. In order 
to calculate accurately the percentage of merged hits, the 
origin of the associated track was shifted to the vertex of 
the event which the trigger particle originated from. 

Previous analyses of low momentum tracks have shown 
that eliminating pairs from both data and mixed events 
with a fraction of merged hits greater than 10% was suf- 
ficient to correct for merging [43|. By discarding pairs 
with more than 10% shared hits we insure that the per- 
centage of merged track pairs is the same in the data and 
the mixed events. 

The correlation function for a given hclicity combina- 
tion of trigger and associated particles was corrected by 
mixed events. After this correction, a small residual dip 
remains, mostly due to track crossing. While the mixed 
events correct for the dip due to true track merging well, 
they do not correct for track crossing as well. The re- 
maining dip is then corrected for using the symmetry 
of the correlations. Since the data should be symmetric 
about A(/) = 0, the data on the same side as the dip are 
discarded and replaced by the data on the side without 
the dip. Then the data are reflected about Ar] = and 
added to the unreflected data to minimize statistical fluc- 
tuations. This method is only applied to \A(f)\ < 1.05 and 
IA77I < 0.67, the region shown in Figure [51 because it is 
computationally intensive and track merging and track 



6 




Act) M Ari 

FIG. 2: (Color online) Two-particle di-hadron correlation function in (Acf), Arj) for 3 < p^'^^^"^ < 6 GeV/c and 1.5 GeV/c 
< p5f=°"=»t<=d ^ pj^rigger Q_i2% central Au+Au collisions at ^/s^ = 200 GeV, not corrected for track merging (left). The 
projection of the correlation function in Acf} for IA77I < 0.042 is shown in the middle and the projection in Ar; for |A0| < 0.17 
on the right. 
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381 
382 
383 
384 
385 



crossing only affect small angular separations. For large387 
Acf) and Arj, the method described in the next section is388 
applied. The track merging correction is done for each389 
di-hadron correlation function separately with the appro- 
priate cuts on p^'*'^"'^, passociatod ^^^^j coUision centrality. 
An example of the di-hadron correlation function beforesw 
and after the track merging correction is shown in Fig- 




FIG. 3: (Color online) The dip region in {Acfi, Arj) uncor- 
rected di-hadron correlations in 0-12% central Au-|-Au col- 

< 6 GeV/c 



trigger 
FT 



lisions at ^/snn = 200 GeV for 3 < 
and 1.5 GeV/c < p--°"-t°d ^ ^trigger ^^^^ 
licity combinations of trigger and associated particles: 

(a) (fttrig, /lassoc) = (1,1), (b) (fttrig , /lassoc ) = (1,-1), (c) 

(/itrig,/iassoc) = ("1,1), and (d) (/itrig,/iassoc) = ("l,"!)- Car- 
toons indicate which dips are from track merging and which 
are from track crossing. 



ure m The slight decrease in the correlation function for 
some data points is an artifact of the correction proce- 
dure and reflects the uncertainty in the correction. 



D. Pair acceptance correction 

With the restriction that each track falls within \r]\ < 
1.0, there is a limited acceptance for track pairs. For Arj 
« 0, the geometric acceptance of the TPC for track pairs 
is « 100%, however, near Arj sa 2 the acceptance is close 
to 0%. In azimuth the acceptance is limited by the 12 
TPC sector boundaries, leading to dips in the acceptance 
of track pairs in azimuth. To correct for the geometric 
acceptance, the distribution of tracks as a function of 
rj and was recorded for both trigger and associated 
particles. A random rj and (f> was chosen from each of 
these distributions to reconstruct a random A77 and Acj) 



for each selection of p^'®^""^. 



5.55 



passociatod^ and centrality. 



5^ 5.5 

<1 

5 5.45 

^ 5.4 
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|5.35 
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^ 5.25 
5.2 



•before 
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-1.5 



-1 -0.5 0.5 

Ari 
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FIG. 4: (Color online) The raw correlation in Arj for di- 
hadron correlations for 3 < p^'^^'^' < 6 GeV/c and 1.5 GeV/c 
< p|?=°"''t<=<i < p!^"sgor Q_^2% central Au-^Au collisions for 
\A(I)\ < 0.78 before and after the track merging correction is 
applied. The data have been reflected about Ar] =0. 
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403 This was done for at least four times as many track pairs453 

404 as in the data and was used to calculate the geometrical454 

405 acceptance correction for pairs. 455 



406 E. Subtraction of anisotropic elliptic flow 45s 

407 background 

460 

408 Correlations of particles with the event plane due t046i 

409 anisotropic flow (U2) in heavy- ion collisions are indirectly 462 

410 reflected in di-hadron correlations and have to be sub-463 

411 tracted for studies of the ridge. This background in A0464 

412 over the interval [—a, a] is approximated by 465 



[—a, a] 



6a0 / dA(j) (1 + 2(4"S) cos (3) 



where a is chosen to be 0.78 so that the majority of the 
signal is contained [H, Ej . The level 6a0 of the back- 
ground is determined using the Zero Yield At Minimum 
(ZYAM) method [3 . The level of the background is taken"' 
as the value of the minimum bin. The systematic errors 
due to the choice of the minimum bin rather than either 

478 

of the two neighboring bins are negligible compared to 
the systematic errors due to the magnitude of V2, dis- 
cussed below. 

481 

The ZYAM method is commonly used for di-hadron 
correlations at RHIC, for example [1, [HI 113], and is jus-""^ 
tificd if the near- and away-side peaks are separated by a""^ 
'signal-free' region. At lower transverse momenta {pT < 2"'"' 
GeV/c) and in central Au-|-Au collisions at ^/snn = 200""^ 
GeV the broadening of the away-side correlation peak""'' 
may cause overlap of the near- and away-side peaks and""^ 
consequently makes the ZYAM normalization procedure 
biased. Alternatively, a decomposition of the correlation 
function using a fit function containing the anisotropic488 
flow modulation of the combinatorial background and 
components describing the shape of the correlation peaks^^^g 
could be used as well psi . |46| . In this paper we use thc^^^ 
ZYAM prescription to remain consistent with our carlier^g^ 
measurements of the near-side ridge [l^]. ZYAM used^^^ 
with conservative bounds on V2 will, if anything, under-^^^ 
estimate the ridge yield. 

For all collision systems and energies studied, the un- 
certainty bounds on V2 were determined by comparing 
different methods for the V2 measurement. We assume 
that the error on the V2 of the trigger and associated 
particles is 100% correlated. Event plane measurements 
of flow and two-particle measurements such as the two- 
particle cumulant method are sensitive to non-flow from 
sources such as jets. These methods may overestimate V2 ■ 
Methods such as the four-particle cumulant method are 
less sensitive to contributions from jets, however, these 
methods may over-subtract contributions from event-by- 
evcnt fluctuations in V2- Therefore these methods under- 
estimate the V2 that should be used for the background5oi 



subtraction in di-hadron correlations [471 ]. For each sys- 
tem at least one measurement which may overestimate V2 
and at least one measurement which may underestimate 
V2 is included. V2 and systematic errors on V2 for Au-|-Au 
collisions at y^sjvw = 62.4 GeV are from comparisons of 
the event plane method using the forward TPCs for the 
event plane determination and the four particle cumu- 
lant method psj . The V2 and systematic errors on 1)2 for 
Au-I-Au collisions at ^/SNN = 200 GeV are as described 
in [HI- The upper bound on V2 is from the event plane 
method using the forward TPCs for the determination 
of the event plane, the lower bound comes from the four 
particle cumulant method, and the average of the two is 
the nominal value. V2 for Cu-|-Cu collisions at ^/snn = 
62.4 GeV and 200 GeV is from [H]. The nominal value 
is given by V2 from the event plane method using the 
forward TPCs for the determination of the event plane 
and the upper bound is from the statistical error for both 
= 62.4 GeV and 200 GeV. For Cu-t-Cu coUisfons 
measurements using the four particle cumulant method 
were not possible due to limited statistics. Instead, for 
Cu-I-Cu collisions at ^/SNN = 200 GeV the lower bound 
is determined by the magnitude of the cos(2A(/)) term 
extracted from fits to p+p data, scaling it by (A'coii), and 
subtracting it from V2 determined using the event plane 
method to estimate the maximum contribution from non- 
flow in Cu4-Cu collisions. For Cu+Cu collisions at ^/snn 
= 62.4 GeV the systematic error is assumed to be the 
same in Cu-|-Cu collisions at both energies. 

With both methods for subtracting the ridge contribu- 
tion to the jet-like yield described below, the systematic 
errors due to V2 cancel out in the jet-like yield, assum- 
ing that V2 is independent of Ar] in the TPC acceptance. 
This assumption is based on the measurements of V2 as 
a function of 77 [HO, [Ml . 



F. Yield extraction 

To quantify the strength of the near-side correlation 
it is assumed that it can be decomposed into a jet-like 
component, narrow in both azimuth and pseudorapidity, 
and a ridge component which is independent of A77. This 
approach is consistent with the method in [l^]. For the 
kinematic cuts applied to and passociated^ 

like correlation is contained within the cuts used in this 
analysis, \At]\ < 0.78 and \A(f>\ < 0.78. 

To study the jet-like correlation and the ridge quanti- 
tatively we adopt the notation from [l2| . and introduce 
the projection of the di-hadron correlation function from 
Eq. (ID) onto the Ar/ axis: 



dN 



dAr] 



dA<j> 



d^N 



dA(j)dAr] ' 



(4) 



and similarly on the A(j) axis: 
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517 
518 
519 
520 
521 

522 

523 



dN 



dA(j) 



|Ar,|e[a,t] dA<j>dAi-i' 



(5) 



To determine the jet-like yield of associated charged 
particles two methods are used. The first method is based 
on Atf) projections. Under the assumption that the jet- 
like yield is confined within IA77I < 0.78, subtracting the 
A(f> projections: 



dN 



dA(j) 



dN 



dA(j) 



0.78 dN 



0,0.78 



1.0 dA(j) 



(6) 



0.78,1.78 



removes both the elliptic flow and ridge contributions. 
Since the second projection is calculated in a larger 
A77 window, it has to be scaled by a factor 0.78/1.0, the 
ratio of the Arj width in the region containing the jet-like 
correlation, the ridge, and the background to the width of 
the region containing only the ridge and the background. 
This subtracts both the ridge and V2 simultaneously since 
within errors both are independent of A77 [l2|, [l^ [sO, IMI ■ 



The jet-like yield Yj 
Eq.® 



is then obtained by integrating 



Y 



A0 



0.78 

dA0 (AcP) 
-0.78 d.A(p 



(7) 



515 The second method for jet-like yield determination is 

516 based on the A77 projection at the near-side: 



dNj 
dA?] 



dN 



dArj 



bAri 



(8) 



-0.78,0.78 



as V2 is independent of pseudo-rapidity within the STAR 
acceptance and therefore only leads to a constant offset 
included in bAri- The background level ^aij is determined 
by fitting a constant background ^a?? plus a Gaussian to 
[Ar]). The yield determined from fit is discarded to 

avoid any assumptions about the shape of the peak and 
instead we integrate Eq. ([8]) over At] using bin counting527 
to determine the jet-like yield Y^^^: 528 



Y 



A)) _ 



dA77 — ^ A77 
-0.78 dAri 



(9) 



The ridge yield Yiidgo is determined by first cvaluat 

then, 



ing Eq. ([5]) over the entire Ar; region to get uii.ci±535 



subtracting the modulated elliptic flow background Ba4,^^'^ 

^At;. 537 

538 

^0.78 

>^ridgo = / dA(f> 



and the jet-like contribution Yj 
dN 



-0.78 



dAcj) 



0,1.78 



SA0hO.78,O.78]-r 



Ai) 



(10)= 
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FIG. 5: Color online) Sample correlations in Arj {\A4>\ < 0.78) 
for 3 < < 6 GeV/c and 1.5 GeV/c < p^==°"^t°d ^ 

pJp"*^'^"" for (a) 0-60% Cu+Cu at = 62.4 GeV (b) 0-80% 

Au+Au at ^/sJTn = 62.4 GeV, (c) 0-95% d+Au at ^/sJTn = 
200 GeV, (d) 0-60% Cu+Cu at ^s]^ = 200 GeV, (e) 40-80% 
Au+Au at ^TsJTn = 200 GeV, and (f) 0-12% central Au+Au 
at y^sjviv = 200 GeV. Lines show the Ar; range where the 
jet-like yield is determined. The data are averaged between 
positive and negative Ar; and reflected in the plot. Shaded 
lines in (f) show the systematic errors discussed in IIIIFI 



mixed event method described in IIIIDI to the standard 
event mixing method and to a sample with a restricted 
z vertex position. The largest difference was seen in the 
central Au+Au data at ^/JWn — 200 GeV. To be conser- 
vative this difference is used as the systematic error for all 
the data sets. The resulting systematic errors are listed 
in Table ini This error is also present for the ridge, since 
the ridge is determined by subtracting Yj. Additionally, 
the systematic error on Yj due to the track merging cor- 
rection does not exceed 1%, the maximum size of the 
correction in the kinematic region studied in this paper. 
This correction does not impact i^idgc- 



G. 2D Fits 



525 We determined the systematic error on the Yj due to un-540 In addition to the standard ZYAM procedure, we also 

526 certainty in the acceptance correction by comparing the54i analyzed the distribution of particles in Eq. ([T|) using 
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FIG. 6: (Color online) Sample correlations in Acjr'^'' 
(|Ar?j < 1.78) for 3 < py'^^^^" < 6 GeV/c and 1.5 GeV/c^"" 
< p-=°"="<=d ^ ^trigger Q_gQy^ Cu+Cu at ^/I^ = 62.4="= 

GeV (b) 0-80% Au+Au at ^s]^ = 62.4 GeV, (c) 0-95%'"^ 
d+Au at ^/s]^ = 200 GeV, (d) 0-60% Cu+Cu at ^/s]^ =56' 
200 GeV, (e) 40-80% Au+Au at ^Sivjv = 200 GeV, and (f)56s 
0-12% central Au+Au at = 200 GeV. Solid lines show 

the estimated background using the ZYAM method with the 
range of V2 used for the determination of the systematic er-^^^ 
rors. The data are averaged between positive and negative 
A(j> and reflected in the plot. Vertical dashed lines show the 
A(j> range where the jet-like correlation is determined. "° 



dAcfxlAi] 



A (l + ^214ACOs(nA0)) 



n.= l 

Yj 



2n(TA(/,,jcrA7j,.j 



with first four coefficients Via, V2A , V3A and V4A of a 
Fourier expansion and a term accounting for tlie jet-like 
correlation on the near-side. This approach is motivated 
by the class of models for ridge production through a tri- 
angular initial condition. If non-flow contributions are 
negligible, V2A corresponds to the average of the prod- 
uct of the trigger particle V2 and the associated particle 
V2 and V3A corresponds to the average of the product 
of the trigger particle V3 and the associated particle W3. 
We use Eq. pT|) to fit the data and extract Vi3A/V2A for 
all collision energies and systems. We allow V3A to be 
negative. A narrow roughly Gaussian away-side peak at 
A0 w TT, which could arise from correlations from the 
production of an away-side jet, would have a negative 
contribution to V3A and a negative V3A could indicate 
that flow is not the dominant production mechanism for 
these correlations. Furthermore, V2A is not constrained 
to the experimental values measured for V2 through other 
means. There is no systematic error on V3A/V2A due to 
the efficiency because any uncertainty in the efficiency 
would change the magnitude of the modulations, given 
by A in Eq. pT|) , but not the relative size of those modu- 
lations, V2A and VsA- The uncertainty due to the fit and 
uncertainty in the acceptance correction is determined by 
fixing the parameters within the range given in Table |IT1 
This gives an uncertainty of <4% on Vsa/V2A- 



IV. RESULTS 



A. Sample Correlations 



571 
572 
573 
574 

542 two-dimensional fits of the form: 

575 

576 
577 
578 
579 

TABLE II: Systematic uncertainties in the acceptance correc-5so 
tion. 581 



trigger 

Pt 

(GeV/c) 


associated 

Pt 

(GeV/c) 


sys. error 
yield 


sys. error 
Gaussian width 


2.0-2.5 


>1.5 


<27% 


<10% 


2.5-3.0 


>1.5 


<18% 


<6% 


3.0-6.0 


>1.0 


<16% 


<6% 


3.0-6.0 


>1.5 


<9% 


<6% 



Figure [5] shows fully corrected Arj projections of sam- 
ple correlations on the near side {\A(j)\ < 0.78) before 
background subtraction for d+Au, Cu+Cu and Au+Au 
collisions at energy ^/sNN = 62.4 GeV and ^/snn = 200 
GcV. The trigger particles were selected with transverse 
momentum 3 < p^'®^""^ < 6 GeV/c and the associated 
particles with 1.5 GeV/c < ^associated ^ p^"gg<=^ The 
data show a clear jet-like peak sitting on top of the back- 
ground and the ridge. The level of the background is 
increasing with energy and system size as expected, as 
more bulk particles are produced in the collision. 

Examples of the complementary projections in A(f> be- 
fore background subtraction for (IA77I < 1.78) are shown 
in Figure [HI The elliptic flow modulated background is 
shown as solid curves. The middle curve corresponds 
to background calculated with the nominal value of V2. 
The upper (lower) curve corresponds to the background 
if the upper (lower) bound on V2 is used instead. Since 
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wc have conservatively assumed that the error on the V2 
of the trigger and associated particles is 100% correlated, 
the background occasionally goes above the signal in Fig- 
ureinif) the away-side. However, since we are focusing 
the near-side we prefer this conservative estimate. Note 
that the uncertainty in the size of the elliptic flow modu- 
lated background affects only the ridge yield but not the 
jet-like yield, since the elliptic flow contribution to the 
jet-like yield in Acj) cancels out in Eq.Q and in Ar] is 
included in 6a?7 in Eq.®. 

Sample background subtracted correlation functions 
1^ from Eq.® and from Eq.® on the near side 

for 3 < py"'^®"' < 6 GeV/c and 1.5 GeV/c < p|?«°"atcd 
< '^^'"^ from Figure [5] and Figure IH] are shown in Fig- 
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FIG. 7: (Color online) Background subtracted sample corre- 
lations for 3 < pJr"^'^'"' < 6 GeV/c and 1.5 GeV/c < p^==°"=it<=d 
< py'^^""^ on the near-side for (a) 0-60% Cu+Cu at ^/sJTn — 
62.4 GeV, (b) 0-80% Au+Au at ^/sWn = 62.4 GeV, (c) 0-95%623 
d+Au collisions at ^/sWn = 200 GeV, (d) 0-60% Cu-f-Cu at624 

= 200 GeV, (e) 40-80% Au+Au at ^/s^ = 200 GeV, 
and (f) 0-12% central Au-|-Au at = 200 GeV. The de 

pendence of the jet-like correlation is shown as a function of 
Acj) {\Ari\ < 1.78) in open symbols and of Arj (| A(j!>| < 0.78) in' 
closed symbols. Lines show the and Ar] ranges where the' 
jet-like yield is determined. The data are averaged between^^' 
positive and negative A77 (At/i) and reflected in the plot. Lines''^'' 
in (f) show the systematic errors on discussed in IIII Fl "^ 



0.25 
0.2 

0.15 
0.1 

0.05 



62 200^ . 
^ d+Au 




-on Cu+Cu 




- • ■ Au+Au 




PYTHIA 6.4.10 Tune A 






m ^ ■ n» * ■ ' 


: 



• ^ • • 



10 



<Npart> 



10^ 



FIG. 8: (Color online) Dependence of jet-like yield on (A'part) 
for 3 < p'r^''''"" < 6 GeV/c and 1.5 GeV/c < < 
^trigger Cu+Cu and Au4-Au at ^^s]7n = 62.4 GeV and 
c?-|-Au, Cu-I-Cu and Au+Au at ^/snn = 200 GeV. Compar- 
isons to PYTHIA, (A^part) =2, are shown as lines. The 5% 
systematic error due to the uncertainty on the associated par- 
ticle's efficiency is not shown and systematic errors due to the 
acceptance correction are given in Table HIl The background 
level and 112 values used for the extraction of these yields are 
given in Table [TVl 



ure [71 For the given kinematic selection, the extracted 
jet-like correlation peaks in both Arj and Acj) projections 
look very similar in all studied systems and collision en- 
ergies. The jet-like yields discussed through the rest of 
the paper are obtained from the Arj projection method; 
the A(j) method is only used for determining the width 
of the jet-like correlation in A^. Below, the dependence 
of the near side jet-like yield and Gaussian width of the 
jet-like correlation peak on collision centrality and the 
transverse momentum of the trigger and associated par- 
ticles are studied in detail. 



B. The near-side jet-like component 



The centrality dependence of the jet-like yield for 3 < 
^tiiggor ^ g GeV/c and 1.5 GeV/c < p|f«°"atod < ptriggcr 

is plotted in FigureH] The jet-like yield at -y/ijviv = 62.4 
GeV is lower than that at ^/sWn = 200 GeV by about a 
factor of three, which can be understood as the result of a 
steeply falling jet spectrum folded with the fragmentation 
function. The measured yields are compared to PYTHIA 
simulations shown as a line in Figure |51 For these stud- 
ies PYTHIA version 6.4.10 [H CDF tune A [H, which 
matches the data from the Tevatron at ^/s =1.8 TeV 
and also describes the pion and proton inclusive spectra 
well at RHIC energies HM^, is used. The PYTHIA 
prediction is somewhat above the data, even for d+Au 
collisions. However, considering the fact that the data 
are from heavy-ion collisions and the PYTHIA simula- 
tions are for p+p collisions, good agreement is unantic- 
ipated. For a given number of participating nucleons, 
(A^part), and collision energy, ^/snn, there is no signifi- 
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FIG. 9: (Color online) Dependence of jet-like yield on p^'*^*"^ 
for 0-95% d-fAu, 0-60% Cu-fCu at ./s^ = 62.4 GeV and 
./sWn = 200 GeV, 0-80% Au+Au at ^/sWn = 62.4 GeV, and 
0-12% and 40-80% Au-hAu at = 200 GeV. Compar- 

isons to PYTHIA are shown as lines. The 5% systematic error 
due to the uncertainty on the associated particle's efficiency 
is not shown and systematic errors due to the acceptance cor- 
rection are given in Table HIl 



cant difference between the d-|-Au, Cu-|-Cu and Au-|-Au 
collisions observed, as expected if the jet-like correlation 
were dominantly produced by vacuum fragmentation. 

The dependence of the jet-like yield on for 
1.5 GeV/c < pl^'^o^ated < p^"gso'- is plotted in Figure [HI 
for all studied collision systems and energies. Data 
from Au-|-Au collisions at y^J/vw = 200 GeV are shown 
separately for peripheral (40-80%) and central (0-12% 



70 J anu ceuirai (u-iz /o) 
Au-f-Au collisions. The jct-like yield increases constantly 
with p^'^^'^' for both Cu-fCu and Au-f-Au and for y^ijvwee? 
= 62.4 GcV and ^i]^ 200 GcV. The effect of aem 
steeper jet spectrum at ^snn ~ 62.4 GeV relative toeeg 
200 GeV discussed above is now reflected in the differ-67o 
ence between the jet-like yields at the two energies. Thise?! 
difference in the jet-like yield between the two studied672 
collision energies increases with p^'^**^"^ from about a fac-673 
tor of two at —2.5 GcV/c to a factor of four at^n 

^trigger ^ g g GeV/c. Comparisons to PYTHIA Simula-"" 
tions are shown as lines in Figure IHl It is surprising how""' 
well PYTHIA is able to describe the p^'®^""^ dependence 



of the jet-like yield in A+A collisions. 



agreement is better at larger (p. 



In general, the 
>4GeV/c), 



while at lower p^'^^""^ values PYTHIA predicts a larger 
jct-like yield than observed in the data. No significant 
differences between d-l-Au, Cu-|-Cu and Au-|-Au collisions 
at ^/snn = 200 GeV are observed. This finding is con- 
sistent with the jet-like correlation arising from fragmen- 
tation. 

The spectra of particles associated with the jet-like cor- 
relation and their comparison to PYTHIA simulations 
for 3 < p^'®^""^ < 6 GeV/c are shown in Figure [TOl For 
the same selection, the mean transverse momen- 

tum fraction zt carried by the leading hadron is larger at 
y^SNN = 62.4 GeV than at 200 GeV due to the steeper 
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10: (Color online) Dependence of jet-like yield on 
^associated ^^j. 3 ^ pt^riggcr ^ g GcV/c for (a) 0-60% Cu-f Cu and 
0-80% Au-I-Au collisions at = 62.4 GeV and (b) 0-95% 

d+Au, 0-60% Cu+Cu, 0-12% Au+Au and 40-80% Au+Au 
collisions at ^snn = 200 GeV. Solid lines through the data 
points are fits to the data. Comparisons to PYTHIA are 
shown as dashed lines. The 5% systematic error due to the 
uncertainty on the associated particle's efficiency is not shown 
and systematic errors due to the acceptance correction are 
given in Table [III 



jet spectrum. This is reflected in softer 



associated 



spectra 



at y/sJiN — 62.4 GeV. The inverse slope parameters from 
an exponential fit to these data arc shown in Table IIIII 
There is no difference seen between Cu+Cu and periph- 
eral Au+Au in either the data points or the extracted 
inverse slope parameter. The inverse slope parameter of 
the central Au+Au data at ^/s/vw = 200 GeV is some- 
what lower than the other data at ^/sNN ~ 200 GeV, 
largely because of the larger yield at the lowest passociatod^ 
This also indicates that there is some modification of the 



TABLE HI: Inverse slope parameters in MeV/c of p^^^°"^^<='i 
spectra from fits of data in Figure [TOl The inverse slope pa- 
rameter from a fit to tt~ inclusive spectra in Au+Au collisions 
[13, m above 1.0 GeV/c is 280.9 ± 0.4 MeV/c for 0-10% 
y/s^ = 62.4 GeV and 330.9 ± 0.3 MeV/c for 0-12% y/s^N 



= 200 GeV. Statistical errors only. 





y/sNN = 62.4 GeV 


y/sNN = 200 GeV 


Au+Au 


332 ± 13 


457 ± 4 (40-80%) 
399 ± 4 (0-12%) 


Cu+Cu 


370 ± 9 


443 ± 3 


d+Au 




438 ± 9 


PYTHIA 


417 + 9 


491 + 3 
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FIG. 11: (Color online) Dependence of the widths in Afji and Aj? on p^r^'^^"' for 1.5 GeV/c < 



associated ^ ^trigger associated 

< Pt ' Pt 



for 



< 6 GeV/c, and (iVpart) for 3 < pf'''°' < 

200 GeV, 0-80% Au+Au at ^^sJTn 



GeV/c and 1.5 GeV/c < p5f=°"-t«d < pt^"sgor ^ gg^^ Q_gQy^ 



o ^ tri 

3 < ' 

Cu+Cu at ^/s]7n = 62.4 GeV and ^/IJTn = 200 GeV, 0-80% Au+Au at ^/sKw = 62.4 GeV, and 0-12% and 40-80% Au+Au 
at y^SNN = 200 GeV. Comparisons to PYTHIA are shown as lines. The 5% systematic error due to the uncertainty on the 



acceptance correction is not shown and systematic errors due to the acceptance correction are given in Table HIl 



677 jet-like correlation at low px- While the agreement with7os 

67S PYTHIA is remarkable for a comparison to A+A colli-709 

679 sions, the discrepancies between PYTHIA and the data7io 

680 are larger at lower momenta and lower energy. This is7ii 

681 expected since PYTHIA is tuned better at higher pT and7i2 

682 higher energy. 713 

683 FigurefTTjshows the Gaussian widths of the Acj) and Arj'"-* 

684 projections of the near-side jet-like peak as a function 

685 of pl^'^ociatod^ g^^^ (A^part) aloug with PYTHIA 

686 simulations. In the most central bin in Au+Au collisions'is 

687 at -y/sivjv = 62.4 GeV it was not possible to extract the 

688 width in Arj because of limited statistics combined with7i6 

689 a residual track merging effect. There are no significant7i7 

690 differences between the widths as a function of p^'^^^'^^^^ 

691 and passociatod f^j. different collision systems except for^^^ 

692 central Au+Au collisions at y^SNN = 200 GeV. While^^^ 

693 no dependence on collision system is observed, there is a^^i 

694 clear increase in the Arj width with increasing (iVpart) in^j^ 

695 Au+Au collisions at -y/ijviv — 200 GeV. This indicatcs^^j 

696 that the shape of the jet-like correlation is modified in^^^ 

697 central Au+Au colhsions at ^/s^n = 200 GeV. PYTHIA^^^ 

698 predicts a greater width in A77 at the lowest p^^°"^^°'^-,2e 

699 than seen in c?+Au or Cu+Cu. 727 

700 Overall it can be concluded that the agreement be-72s 

701 twecn the different collision systems and energies shows729 

702 remarkably little dependence of the jet-like per trigger73o 

703 yield on the system size. In contrast to the periph-731 

704 eral Au+Au data, the central Au+Au data show in-732 

705 dications that the jet-like correlation is modified. The733 

706 model in [59j , a hypothesis for the formation of the ridge734 

707 through gluon bremsstrahlung, does not produce a ridge735 



broad enough to agree with the data, however, it is possi- 
ble that a similar mechanism could explain the broaden- 
ing of the jet-like correlation. Similarl y rn odels for ridge 
production by turbulent color fields j27l l60j predict a 
broadening of the jet-like peak in Arj which is not wide 
enough to describe the ridge but may explain the data in 
Figure [TT] 

C. The near-side ridge 

In [l2|, we reported detailed studies of the ridge in 
Au+Au collisions at ^/snn = 200 GeV function 
of p^'^®'''^ and passociatcd^ Here we investigate the ridge 
centrality, energy and system size dependence. The de- 
pendence of the ridge yield on (A^part) for Cu+Cu and 
Au+Au collisions is shown in Figure [12] for both ener- 
gies studied. Table HV] shows 6a0 values and V2 of trig- 
ger and associated particles for all collision systems and 
energies studied in Figure [121 The centrality bins are 
characterized by the fraction of geometric cross section 
a/CTgco, average number of participants (A'part) and num- 
ber of binary collisions {Ncoi\)- Contrary to the jet-like 
yield, which shows little dependence on centrality, the 
ridge yield increases steeply with (A'part)- Within er- 
rors, there is no difference in ridge yield between Cu+Cu 
and Au+Au collisions at the same (A'part) at a given en- 
ergy, demonstrating the system independence of the ridge 
yield. 

The energy dependence of the ridge yield is potentially 
a sensitive test of ridge models. Comparing the two col- 
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System 


C /Cgco 


(Afpart) 


(A^coll) 




^^trigger^ 




V2A 


V3A 




[%] 








[%] 


[%] 


[%'] 


[%'] 


d+Au 


MB 


8.3+0.4 


7.5+0.4 


0.158+0.006 


0.0 


0.0 


69 + 1 


-285 + 14 


200 GeV 


















Cu+Cu 


0-10 


96+3 


162+14 


0.885+ 0.010 




7 0+1.9 
1 "3-0.7 


84 + 5 


54 + 5 


62.4 GeV 


10-30 


64+1 


87.9+-7.9 


0.515+ 0.006 






228 + 5 


5 + 5 




30-60 


25.7+0.6 


27.6+1.6 


0.201+ 0.005 




11.8l«;? 


471 + 10 


-77 + 10 


Au+Au 


0-10 


320+5 


800 + 74 


3.582+0.019 


8.5t-j 




63 + 3 


42 + 2 


62.4 GeV 


10-40 


169+9 


345 + 44 


1.846+0.010 


i7.5i«:^ 




258 + 2 


57 + 2 




40-80 


42+8 


51 + 16 


0.446+0.009 




18.6l3;4 


456 + 8 


-33 + 8 


Cu+Cu 


0-10 


99+1 


189+15 


1.759+ 0.007 




».»_2 2 


128 + 2 


65 + 1 


200 GeV 


10-20 


75+1 


123.6+9.4 


1.206+ 0.006 




11.31?;° 


212 + 2 


62 + 2 




20-30 


54+1 


77.6+5.4 


0.815+ 0.006 






296 + 3 


47 + 3 




30-40 


38+1 


47.7+2.8 


0.540+ 0.006 


10.71^6 


13.01";^ 


380 + 4 


9 + 4 




40-50 


26.2+0.5 


29.2+1.6 


0.337+ 0.005 




12.6t»- 


506 + 6 


-37 + 6 




50-60 


17.2+0.4 


16.8+0.9 


0.209+ 0.005 


10.3l°;« 


■19 o+O.O 
lz.o_2 g 


657 + 10 


-162 + 10 


Au+Au 


0-12 


316+6 


900+71 


13.255+0.003 




7 q+1.5 


81 + 1 


49.9 + 0.2 


200 GeV 


10-20 


229+5 


511+34 


4.683+0.007 


14.6i-° 


is.otli 


225 + 1 


67.1 + 0.6 




20-30 


164+5 


325+23 


3.222+0.006 


i7.8t^:i 




344 + 1 


71.2 + 0.8 




30-40 


114+5 


199+16 


2.094+0.006 


18.9l^i 


is.ilil 


430 + 1 


68 + 1 




40-50 


76+5 


115+12 


1.267+0.006 




19.01^:^ 


461 + 2 


57 + 2 




50-60 


48+5 


61+8 


0.738+0.006 


lo.il'J 


17.61^:^3 


478 + 3 


18 + 3 




60-70 


28+4 


30+5 


0.386+0.006 






549 + 6 


-25 + 6 




70-80 


15+2 


14+3 


0.183+0.006 






754 + 12 


-195 + 12 



TABLE IV: The background terms bA4> (see Equation IIII E|) . elliptic flow values of trigger (dj'^'^^'^') and associated particles 
^^assoc^j^ and Fourier coefficients V2A and Vsa from 2D fits for different collision energies, systems and centrality bins defined 
by the fraction of geometric cross section (a/CTgco), average number of participants ((Ai'part)) and binary collisions ((A'^coii)) for 
the data in Figure [S] 



736 lision energies studied, the ridge yield is observed to bc752 

737 smaller at ^/snn = 62.4 GeV than at ^/snn = 200 GeV. 753 

738 Similar behavior was also observed for the jet-like yield. 754 

739 Therefore a closer investigation of the centrality depen-755 

740 dence of the ratio yridgc/Xjot is reported in Figure 1131 756 

741 The ratio of the yields is independent of collision energy757 

742 within errors. For the same kinematic selections, the data758 

743 at y^SNN = 62.4 GeV correspond to a lower jet energy759 

744 which may imply the decrease of the ridge yield with the76o 

745 parton energy, as observed for the jet-like yield. 76i 

762 

746 A recent STAR study of the ridge using two-particle az-^^^ 

747 imuthal correlations with respect to the event plane [l^ 

748 in Au+Au collisions at ^/snn = 200 GeV shows that the764 

749 ridge yield is dependent on the angle of the trigger par-765 

750 tide relative to the event plane. Another STAR study766 

751 of three-particle correlations in pscudorapidity [l3| in767 



Au+Au collisions at y'J/vAr = 200 GeV shows that within 
current experimental precision particles in the ridge are 
uncorrelated in A77 with both the trigger particle and 
each other. These observations along with the PHOBOS 
measurement of the ridge out to Ar] =4 [l^ provide sub- 
stantial experimental constraints on theories for the pro- 
duction of the ridge. In addition, the CMS experiment 
recently observed the ridge in high multiplicity p+p col- 
lisions at ^/s =7 TeV [61|. We consider these results in 
addition to the results presented in this paper in order 
to evaluate our current theoretical understanding of the 
ridge. 

The momentum kick model (l9| with the same kine- 
matic selection criteria applied to charged particles de- 
scribes the increase of the ridge yield with centrality 
quantitatively in Au+Au collisions at ^/sWn = 200 GeV. 
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768 For a given collision energy the ridge yield for Cu+Cu 

769 collisions is predicted to approximately follow that for 

770 Au+Au collisions for the same number of participants, 

771 however the prediction for Cu+Cu collisions is systemat- 

772 ically above that for Au+Au collisions at both energies. 

773 This difference between the two systems is not corrob- 

774 orated by the data. For the same nucleus-nucleus col- 

775 lisions at different energies, the ridge yield in the mo- 

776 mentum kick model is predicted to scale approximately 

777 with the number of medium partons produced per partic- 

778 ipant which increases with increasing collision energy as 

779 (luy/s)^ [13 ■ According to this prediction, Yridgo should 

780 increase by a factor of 1.6 in ^/snn = 200 GeV rela- 

781 five to 62 GeV. This is in agreement with our measure- 

782 ment. We note that the momentum kick model is likely 

783 to have difficulty describing the observed dependence of 

784 the ridge yield on the event plane, since it would likely 

785 predict larger ridge out-of-plane. It would also likely have 

786 difficulty explaining the absence of correlations between 

787 particles in the ridge and the jet-like correlation. 

788 The model where the ridge arises from the coup ling 

789 of induced giuon radiation to the longitudinal flow [18l|802 

790 is in qualitative agreement with the observed increasesos 

791 of l^idgc as a function of (A'part) since the size of the804 

792 ridge should depend roughly on the average path lengthso5 

793 traveled by a hard parton. However, it is not obvioussoe 

794 that this model can describe the collision system andgo^ 

795 energy dependence of ridge yield reported in this paper. 

796 Moreover, the large extent of the ridge in Ar] and absenceg„g 

797 of correlation among ridge particles clearly disfavors thiSgj„ 

798 physics mechanism for ridge formation. 

799 The radial flow plus trigger bias model [lH, [2^ ^^su 

800 would predict an increasing ridge yield with increasingsu 

801 A^part- This model would predict a larger ridge in-planesw 
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FIG. 12: (Color online) Dependence of the ridge yield on^^^ 
(A^part) for 3 < p'r^'"""' < 6 GeV/c and 1.5 GeV/c < p5f=°"="<=d33, 
< py"'^'^'^'' for Cu+Cu, and Au+Au at ^/s^Tn = 62.4 GeV and833 
Cu+Cu and Au+Au at ^/snn ~ 200 GeV. Systematic errors834 
due to V2 are shown as solid lines. The 5% systematic error835 
due to the uncertainty on the associated particle's efficiency 
is not shown and systematic errors due to the acceptance^^j 
correction are given in Table HIl 




FIG. 13: (Color online) Ratio of the ridge and jet-like yields as 
a function of (A^part) for 3 < py"*^*^"" < 6 GeV/c and 1.5 GeV/c 
< p^ssociatcd ^ pl^riggcr gygtematic errors due to V2 are shown 
as solid lines and systematic errors due to the acceptance 
correction are given in Table HIl 



because of the larger surface area in-plane, in agreement 
with the data. Since the ridge in this model arises from 
medium partons, particles in the ridge are not expected 
to be correlated with each other, again in agreement with 
the 3-particle A77 correlation data. 

The mechanisms for the production of the ridge 
through glasma initial state effects (28l430j are able to ex- 
plain the observed large A77 extent of the ridge. However, 
it is not obvious what this class of models would predict 
for the collision system and energy dependence discussed 
in this paper as well as for other ridge properties includ- 
ing its p^'^^"'^ and ^associated dependence, the event plane 
dependence, and the absence of correlated structures in 
three-particle A77 correlations. While there are calcula- 
tions of untriggered di-hadron correlations for the glasma 
model, there are no calculations to compare to high-pT 
triggered correlations. It should be noted that these cal- 
culations not only include the glasma initial state but 
also hydrodynamical flow so some of the ridge in these 
models is created by flow. If the ridge is produced by 
the same mechanism in p+p and A+A collisions, models 
where the ridge is produced by initial state effects such 
as the glasma model may be the only models able to ex- 
plain both the p+p and A+A data simultaneously since 
hydrodynamical flow in p+p collisions is expected to be 
small if not negligible. Therefore quantitative calcula- 
tions of the ridge in this class of models and the identified 
particle spectra measurements in high multiplicity p+p 
collisions are essential for understanding the production 
mechanism of the ridge. 

Models describing the ridge in terms of quadrupole, tri- 
angular and higher order components, from initial ec- 
centricity fluctuations predict a ridge yield that increases 
with A^pait m, m, [ssl. ItoI in qualitative agreement with 
data. Motivated by these models, we applied two dimen- 
sional fits to our data with a two dimensional Gaussian 
to describe the jet-like component and Arj independent 
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FIG. 14: (Color online) 1/3a/V'2A ratio as a function of (A^part) for 3 < < 6 GeV/c and 1.5 GeV/c < p^==°"=>t<=d ^ pj^^iggor 

Statistical errors only. The systematic error of <4% not shown. The data are compared to hydrodynamical calculations using 

(a) lattice equation of state with three different freeze-out temperatures in Au+Au collisions at y^ijvjv = 200 GeV from [6^ . |65| ] 

(b) using a 3+lD viscous hydrodynamical model with three different viscosity to entropy ratios, 77/s [66l. l67|. and (c) a 3+ ID 
model incorporating both hydrodynamics and transport model |6^. (a) and (c) use the same kinematic cuts as the data 
while (b) uses 1.5 < pr < 4.0 GeV/c for both p^''^^'^'^'' and p^===°<=i'>-t«d 
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terms given by Eq. ffTT]l. Figure fT4l shows V3A/V2AB75 
as a function of (A^part) from these fits. The values of876 
the Fourier coefficients V3A £^nd V2A are given in Ta-877 
blellVi We allow V3A to be negative and V3A is negative878 
for (i+Au collisions and peripheral A-\-A collisions. An879 
approximately Gaussian peak from an away-side jet-likcsao 
correlation would give a negative V3A ^^nd indicate thatsai 
V3A is dominantly from non-flow. The ratio V3a/^2A882 
evolves from negative values in d+Au and peripheral883 
A-\-A collisions to positive values at larger (A^part)- 884 
Positive values for V3A are consistent with expect at ionsB85 
from triangular initial conditions. Contributions to a886 
jet-like peak on the away-side, such as that observed in^s? 
d-hAu, would lead to V3A/V2A < {vzjv^f' because V3A888 
would be an underestimate of v\ and V2A would be anaag 
overestimate of v\. For both Au-hAu and Cu-hCu colli-89o 
sions, the ratio is independent of collision energy withiuggi 
errors and is largest in the most central collisions with892 
significant deviation between the two colliding systems ats93 
the same (A^part)- It is not clea.r at this point whether thes94 
representation of the data in Figure [T^ where the vi hasggs 
been subtracted, or in Figure [HI where the V2 is explic-Bse 
itly included, gives the most insight into the productiongg? 
mechanism of the ridge. sss 
We compare the data to three hydrodynamical modelssw 
in Figure [T31 noting that these data are in a momen-soo 
tum range approaching the limit where hydrodynamicalsoi 
models are expected to be vaUd. By comparing differcnt902 
models to the data we are able to see whether the data^oa 
can constrain the initial state. In Figure [Hk we comparc<>04 
the data to 2+ ID hydrodynamical model calculations for^os 
Au+Au collisions at ^/snn = 200 GeV for various hadronsoe 
resonance gas freezeout temperatures and using a mod- 



ified Glauber initial state (6J|. These calculations use agoa 



lattice equation of state and a viscosity to entropy ratiogog 
rj/s = l/47r [g^. We note that this model does not in-gio 



elude resonance decays. These predictions agree with the 
data for all but the three most peripheral centrality bins 
for a hadron gas freezeout temperature of 170 MeV. 

In Figure [Mb we compare the data to 3+lD hydrody- 
namical model calculations for several values of 7?/s and a 
kinetic freezeout temperature of 130 MeV [6^, [g^]. This 
model uses a Glauber initial state modified to generate 
structures like those expected from flux tubes and dif- 
fers from [g^l because it includes nontrivial longitudinal 
dynamics and the effect of varying the ry/s at a given 
freeze-out temperature. The calculations in (66l.[67| were 
limited in momentum so it was only possible to compare 
to calculations with momenta of 1.5 < pr < 4.0 GeV/c 
for both py'ss" and p^ssociatod^ Quahtatively this model 
reproduces the trends observed in the data and it tends 
to support a higher value of 77/s. 

Figure [Tib shows a comparison of the data to a 3+ ID 
event-by-event transport -I- hydrodynamical model cal- 
culation (gsI . [59} . In this model, the initial conditions 
such as long range rapidity correlations and fiuctuations 
in the transverse energy density profile are provided by 
the UrQMD model [Hi ^ . The hydrodynamic evolu- 
tion starts at 0.5 fm. A transition from hydrodynamic 
evolution to the transport approach is followed by final 
state rescatterings and resonance decays. The predic- 
tions shown in the figure are scaled by a factor 0.5. Re- 
quiring a high pt trigger particle skews the calculation 
towards events with a hot spot in the initial density pro- 
file, which may lead to a preference for events with a high 
^3 value. This model qualitatively describes the (iVpart) 
dependence of V3A/V2A for Au-|-Au collisions at both 
y/IWN = 62.4 GeV and 200 GeV. 

Despite differences in initial conditions, transport pa- 
rameters and freeze-out requirements, all three models 
shown are able to reproduce the qualitative trends of 
V3A/V2A versus (iVpart) for large (iVpart)- Since V3A is 
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911 expected to be largely independent of centrality [s^, iTSjgss 

912 this likely reflects the model's accuracy in predicting cen-^^^ 

913 trality dependence of V2A • The models in I65L 16611 repro-g55 

914 duce the data fairly well quantitatively, but |65| achieves„„ 

915 the best fit by varying the temperature while [66[ varieSgj^ 

916 the viscosity. Since they also have different initial condi-gjg 

917 tions, wc infer that agreement with the data is possible^jg 

918 even with different assumptions and parameters duringg^^ 

919 the hydrodynamical evolution. It will be interesting tOg^^ 

920 see if such similarity persists at lower transverse momentagj^^ 

921 where the hydrodynamic calculations are more reliable. g^j 

922 These studies imply that while fluctuations of the initialgg^ 

923 state are needed to induce odd higher order w„ terms, g^s 

924 the observable V3A/V2A is rather insensitive to the exactg^u 

925 details of the model. This measurement does howeverg^, 

926 enforce added restrictions to model implementations andg^g 

927 should therefore be added to the suite of results, suchg^g 

928 as identified particle spectra, yields and px dependencCg^^ 

929 of U2, currently used to validate theories. Future the-g^^ 

930 oretical and experimental studies will be needed to de-g^j 

931 tcrminc whether these models are sufficient to describCg^j 

932 the complete ridge, and/or whether there are substantial 

933 non-flow contributions to the Fourier coefficient V3A in 

934 central A+A collisions. 



by vacuum fragmentation. 

The ridge is observed not only in Au+Au collisions 
at y^J/vlv = 200 GeV as we reported earlier, but also in 
Cu+Cu collisions and in both studied collision systems 
at lower energy of ^/snn = 62.4 GeV. This demonstrates 
that the ridge is not a feature unique to Au+Au collisions 
at the top RHIC energy. We observe two trends which 
set significant limits to models. First, when the ridge 
is measured using the standard ZYAM model, the ridge 
is comparable in Cu+Cu and Au+Au collisions and the 
energy dependence of the ridge is the same as the energy 
dependence of the jet-like correlation. Second, when we 
subtract the jet-like correlation and calculate the third 
component of the Fourier decomposition V3A, which is 
v| in the absence of non-flow, we see different trends for 
Cu+Cu and Au+Au but no difference between the two 
energies. The combination of these data with future mea- 
surements at lower RHIC energies {y/sNN = 7-39 GeV) 
as well as studies at the LHC will therefore be a pow- 
erful tool for the distinction between various theoretical 
models for the production of the ridge. 



935 V. CONCLUSIONS 

974 

936 The energy and system size dependence of near-side 

937 di-hadron correlations enables studies of the jet-like cor-975 

938 relation and the ridge at flxcd densities with differcnt976 

939 geometry. The reasonable agreement of the jet-like cor-977 

940 relation with PYTHIA is surprising, especially consider-978 

941 ing that PYTHIA is a p+p event generator. There is979 

942 remarkably little dependence on the collision system at98o 

943 both ^/SNN = 62.4 GeV and 200 GeV except for cen-98i 

944 tral Au+Au collisions at y^SNN = 200 GcV. In central 

945 Au+Au collisions at ^/sNN = 200 GeV the jet-like cor-983 

946 relation is substantially broader and the spectrum softer984 

947 than in peripheral collisions and than those in collisions985 

948 of other systems in this kinematic regime. This may indi-986 

949 cate that fragmentation is modiflcd in these collisions SO987 

950 that the parton fragments softer, perhaps due to a mech-9a8 

951 anism such as gluon bremsstrahlung. This indicates that9a9 

952 the near-side jet-like correlation is dominantly produced99o 
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